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ABSTRACT 

We submit that non thermalized support for the outer intracluster medium in relaxed galaxy clusters is provided by turbulence, driven 
by inflows of intergalactic gas across the virial accretion shocks. We expect this component to increase briskly during the cluster 
development for z < 1/2, due to three factors. First, the accretion rates of gas and dark matter subside, when they feed on the outer 
wings of the initial perturbations in the accelerating Universe. Second, the infall speeds decrease across the progressively shallower 
gravitational potential at the shock position. Third, the shocks eventually weaken, and leave less thermal energy to feed the intracluster 
entropy, but relatively more bulk energy to drive turbulence into the outskirts. The overall outcome from these factors is physically 
modeled and analytically computed; thus we ascertain how these concur in setting the equilibrium of the outer intracluster medium, 
and predict how the observables in X rays and ^/waves are affected, so as to probe the development of outer turbulence over wide 
cluster samples. By the same token, we quantify the resulting negative bias to be expected in the total mass evaluated from X-ray 
measurements. 

Key words, galaxies: clusters: intracluster medium — turbulence — X-rays: galaxies: clusters — methods: analytical. 



1. Introduction 

In galaxy clusters the gravitational potential wells set by Dark 
Matter (DM) masses M ~ 1O 15 M are filled out to the virial 
radius R ~ Mpc by a hot thin medium at temperatures k%T ~ 
several keVs, with central particle densities n ~ 10~ 3 cm -3 . 

Such a medium constitutes an optimal electron-proton 
plasma (that we appropriately name IntraCluster Plasma, ICP), 
with its huge ratio of the thermal to the mean electrostatic en- 
ergy ksT le 2 n 1 ^ ~ 10 12 , and the relatedly large number nA i D = 



(k B T/47re 2 n l/3 ) 3/2 ~ 10 16 of particles in the Debye cube. 



It 



emits copious X-ray powers Lx oc n L T l/Z ~ 10 45 erg s 1 via 
thermal bremsstrahlung, but with long radiative cooling times 
over most of the cluster volume. 

Thus on scales longer than the electro-proton mean free path 
the ICP constitutes a quasi-neutral, simple fluid with 3 degrees of 
■ freedom in thermal equilibrium, and with effective particle mass 
. £h ' ^ m p ~ 0-6 m p m terms of the proton's m p . This affords precision 
, modeling on the radial scale r - to begin with - for the dis- 
j_i ■ tributions of density n(r) and temperature T(r), as to match the 
C$ j wealth of current and upco ming data concerning the emissions 
in X r ays (e.g., reviews bv lSnowden et al.ll2008t iGiacconi et al.l 
2009), and the strengths y oc n T e of the Sunyaev-Zel'dovich 
(1972) scattering in /iwaves (e. g., lBirkinshaw & Lanca ster 2007; 
ISchafer and Bartelmannll2007l) . 

In fact, simple yet precise modeli ng is provided by the 
Supermodel (SM; ICavaliere et al1l2009l) . This is based on the 
run of the ICP specific 'entropy' (adiabat) k = ksT/n 2 ^ set by 
the processes for its production. 

The entropy is raised at the cluster centers due to the en- 
ergy discharged by AGN outburs t s (see iVaj ageas & Silk 1999; 
Wu et al. 200C| ICavaliere etafl 120021: iMcNamara & N ulsen 
2007b. and by deep merge rs (see iMcCarthv et al.l 120071; 



ing (see IZuHone et al.l 12010). At the other end, much entropy 
is continuously produced at the virial boundary; there the ICP is 
shocked by the supersonic gravitational inflo w of gas accreted 
from the environment along with the DM (see lTozzi & Normanl 
|200TI lvoi^l2005l;lLapi et al. 2005), and is adiabatically stratified 
into the DM potential well. 

These physical processes concur to originate ICP entropy 
distributions with spherically averaged profiles k(r) - k c + (kR - 
k c )(r/R) a ; these comprise a central floor k c « 10 - 100 keV 
cm 2 , and an outer ramp with slope around a « 1, rising to ad- 
join the boundary values kR ~ some 10 3 keV cm 2 . Such values 
and shapes are consiste nt w ith recent o bservational analyses by 
ICavagnolo etail d2009l) and lPratt et aLld2010l) out to r « R/$\. 

The ensuing gradient of the thermal pressure p(r) oc 
k(r) n 5 ^ 3 (r) is used in the SM to balance the DM gravitational 
pull -G M(< r)/r 2 and sustain hydrostatic equilibrium (HE) 
out to the virial boundary. The HE equation may be written 
as dlnT/dlnr = 3 a/5 - 2b/5 in terms of the entropy slope 
a(r) = dln£/dlnr, and of the potential to thermal energy ra- 
tio b(r) = fmi p v 2 /kBT with v 2 = G M(< r)/r. Whence we di- 
rectly derive the temperature profile (see ICavaliere et al.l [2009. 
their Eq. 7) 
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in terms of the entropy run k(r) and the boundary values at r = R. 
The density follows n(r) = [£Br(r)/&(r)] 3 '' 2 , so that T(r) and 
n(r) are linked, rather than independently rendered with multi- 
parametric expressions as in other approaches. From T(r) and 



Markevitch & Vikhlininl 120071) often followed by inner slosh- 



1 Note that R/2 w ^joo ~ 2i?20o/3 holds in terms of the radii in- 
side which the average DM overdensity relative to the critical universe 
amounts to 500 and 200, respectively. 
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n(r) the X-ray and SZ observables are readily derived and com- 
pared with data. 

In preparation to developments given below, we stress that 
the few parameters specifying k(r) are enough for the SM to 
provide remarkably good fits to the detailed X-ray data on sur- 
face brightness and on temperature profiles of many clusters (see 
iFusco-Fem iano et al] |2009l) . These include central temperature 
profile s of both main classes identified by Mole ndi & Pizzolatol 
(120011) : the cool-cored CCs with a central dip, and the centrally 
flat, non-cool-cored NCCs. The SM intrinsically links these mor- 
phologies to low or high entropy levels k c ~ 10 1 or ~ 10 2 
keV cm 2 , respectively; these are conceivably imprinted by en- 
ergy inputs from central AGN outbursts or from deep mergers 
dCavaliere et al.ll2009llFusco-Femiano et alj|2009t) . 

The SM also covers diverse outer behaviors including 
cases where t he entropy pro duction decreases and its slope 
a abates (see lLapi et all 120101 and data references therein), 
to the effect of producing steep temperature profiles. The in- 
terested reader may try for her/himself other clusters on us- 
ing the fast SM algorithm made available at the website 
http://people.sissa.it/~lapi/Supermodel/. Here we pursue another 
consequence of the diminishing entropy production, namely, tur- 
bulence arising in the outskirts of relaxed CC clusters. 

A num ber of observations (in particular with the Suzaku 
satellite, see iGeoree et alJ20 09. Bautz et al.l2009l l Hoshino et al. 
|20 1 Ol iKaw aharada et al. 201 Op, support the notion that HE may 
be contributed by non-thermalized, turbulent motions occurring 
on scales of several 10 2 kpc inwards of R. On the other hand, sev- 
eral simu lations resolve a variety of shocks in and around clus- 
ters (e.gjRicker & Sarazinl200l[|Markevitch & Vikhlininl2007l 
ISkillman et al.ll2008ilVazza et alJl2010h . 

We focus on the accretion shocks that originate at the virial 
boundary from inflows of gas preheated at temperatures 10 6 K 
by sources like outer stars and AGNs or by hydrodynamical pro- 
cesses like shocks around filaments. The above two pieces of in- 
formation lead us to investigate whether the physics of the virial 
accretion shocks indeed requires turbulence to develop also in 
the outskirts of relaxed clusters under the smooth inflows that 
preva il there dFakhouri et al. 2010, lGenel et al.l2010llWang et al.l 
2010). 



2. Virial accretion shocks: entropy and turbulence 

The key physical agent is constituted by the residual bulk flows 
downstream the virial accretion shocks. The latter actual ly form 
a complex network ('shock layer', see La pi et al]|2005l) modu- 
lated to different strengths by the filamentary structure of their 
environment. Shock curvature leading to baroclinic instabilities 
and vortical flows, and/or sheared inflows inevita bly arise as pi- 
oneere d in the context o f cosmic structures by iDoroshkevichl 
(1 19731) and lBinnevI (1 19741). and recently demonstrated by several 
hydro-simulations (see llapichino & Niemeverll2008i iRvu et al~l 
l20MlI^uetaT]l2009ll^detal.ll2010l) . " 
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2. 1 . Physics of virial accretion shocks 

Boundary shocks of high average strength arise in conditions 
of intense inflows of outer gas with supersonic speed v\ corre- 
sponding to Mach numbers M 2 = v 2 /c 2 ~ 10 2 relative to the 
outer sound speed c\ = (5 ksT[/3 /-im p ) 1 ^ 2 . These shocks effec- 
tively thermalize the inflows, to produce postshock temperatures 
close to the ceiling ks Ti ~ ymi p v 2 /3 that would mean full con- 
version of the infall m p v 2 /2 to thermalized energy 3 kg T/2fi per ^4> - 



Fig. 1. To illustrate the gravitational potential governing gas in- 
fall. As the DM outskirts develop (see Sect. 2), both the virial R 
and the turnaround radius R ta shift outwards, retaining the ratio 
Rtd/R ~ 2. Meanwhile, the outer potential becomes shallower 
and AO lower from the value marked in black to that in blue. 
Moreover, the shock position R s slowly outgrows R, lowering 
yet the drop to the value marked in red. 



electron-proton pair. But even strong shocks hovering at R s « R 
lea ve some residual postshock bulk flows with speed \>2 ~ vi/4, 
see lLapi et al.1 (120051) . These correspond to a kinetic energy ratio 
v 2 /v 2 *6.3%. 

Such a ratio is bound to grow, however, during the 
outskirts development. On the DM side, the latter develop 
insid e-out by secular accretion after the early central collapse 
(see lLapi & Cavaliere 2009, Wang et aTl l2010i and references 
therein). Such a trailing accretion feeds scantily on the outer, de- 
clining wings of the initial DM density perturbation that devel- 
ops into a cluster, and is further impaired b y the cosmic expan - 
sion accelerated by the Dark Energy (cf. iKomatsu et al.l l2010). 
In these conditions, the inflows will peter out, shock thermaliza- 
tion will be reduced and eventually the shock themselves weak- 
ened, to leave postshock bulk energies enhanced well above the 
ratio 6.3%. 

To quantify the issue, we describe the perturbation shape in 
terms of the effective powerlaw 6M/M oc M~ e that modulates the 
mass excess SM accreting onto the current mass M; in particular, 
low values of the shape parameter e < 1 apply to t he perturbatio n 
body, but e grows larger for the outskirts (see iLu et al.l 12006). 
A shell 6M will collapse on top of M when 6M/M attains the 
critical threshold 1.69 D~ l (t) in terms of the linear growth factor 
D(t). So the parameter e also modulates the average mass growth 
reading M(t) oc D 1 ^ oc tr^ e , with the growth factor represented 
by the powerlaw D(t) oc f 1 in terms of the exponent d; the latter 
decreases from 2/ 3 to 1/2 as the r edshift lowers from values 
z ^ 1 to z < 1/2, cf . I Weinberg! d2008l) . Thus the outskirts develop 
at accretion rates M/M ~ d/et that lower as e takes on values 
exceeding 1 in the perturbation wings (and formally diverging in 
voids), and as d decreases to 1 /2 at late cosmic times. 

On the ICP side, the outer gas will accrete with a lower infall 
speed vi. As illustrated in Fig. 1, the latter is set to v\ — 2 AO by 

the outer gravitational potential drop AO = J R drG 5Mj r ; this 
is experienced by successive shells of DM and gas that - after an 
initial expansion - turn around at the radius R ta »2fi to begin 
their infall toward the shock at R s ss R. Thus the potential drop 
A(f> = AO/ v\ (normalized to the circular velocity vi = G M/R at 
r - R) reads 



1 - Cg/gg) 
3e-2 



3e-2 



(2) 



A. Cavaliere et al.: Non Thermal Support for the ICM 



3 



Non-thermolized support 




1 .0 



Fig. 2. Non-thermalized pressure support 6{r) = p n th//>th nor- 
malized to the boundary value 6r computed after Eq. (8). The 
blue line refers to i — 0.9 R, red line to t — 0.5 R, and green line 
to{ = 0.2/?. 



This is seen (cf. Fig. 1) to become shallower during the outskirts 
development as e exceeds 1 ; then the approximation A<p « (3 e - 
2)- 1 * (3e) _I applies. 

Actually, the shock position R s slowly outgrows the virial R 
to approa ch j? t a, and this yields an even lowe r effective potential 
drop (see IVoitet al.l 12003b lLapi etalj|2005h . In fact, it can be 
shown that for e > 1 the shock posit i on ma y be approximated as 
R s ss 2R(l -4e~ 2 ), see lLapietaLld2010l) ; so from Eq. (2) we 
obtain * 4 e~ 2 . 

The direct and indirect effects of M dwindling combine into 
the accretion rate scaling M oc v\M/R oc v 3 /A(p; whence the 



infall speed follows 



vi cc M l/3 (A(f>) 



,1/3 



(3) 



This is indeed reduced strongly during the late development of 
the outskirts when both M oc d/e subsides and A<p oc e~ 2 lowers, 
so that the overall scaling vi oc d 1 ^ 3 /e applies for e > 1; specif- 
ically, as € increases from 1 to 1.5 and then to 3, the prevailing 
Mach numbers decline from Ai 2 ~ 10 to 6 and then to 3. 

Such values are consistent with the Mach number distribu- 
tions at low z sliced for flows of preheated g as into the clus- 
ter, as found by num erical simulations (e.g.. iRvu et al. 1 12003: 
Skillman et al. 2008, see their Fig. 6 and Sect. 4; see also 
Vazza et alj|2010h . 



2.2. Weakening shocks and entropy demise 

With vi lowering toward transonic values the shoc k strength 
will e ventually weaken. We recall (see Appendix B in lLapi et alj 
120051) that the Rankine-Hugoniot conservation conditions for a 
standing shock give the postshock temperatures and densities in 
the form 



Ti 8 + 16 16 M 2 ' «i 3 + M 2 



(4) 



As the Mach number decreases from values Ai > 3 (defining 
strong shocks) to transonic M 2 ~ 1, the temperature ranges 
widely from ksTo - 3 fim p v 2 /16 to T2 — T\, while the den- 
sity varies mildly from «2 — 4ni to «2 — «i- Correspondingly, 



the shock-generated entropy k^ = ksT2/r£ drops from values 
5 M 2 k\/4 2 ^ 16 ~ some 10 3 keV cm 2 typical of strong shocks to 
intergalactic values k\ = ksT\/n 2 ^ 3 < 10 2 keV cm 2 . 

By the same token, the entropy out er ramp will abate. Its 
slope at the boundary has been derived by Cavali ere et al.l {2009) 
from the jumps at the shock and the adjoining HE recalled in 
Sect. 1, to read 



a R = — (5 - b R ) ; 



(5) 



here ft. = 4(1 + e/d)/[5 + 2(1 + e/d)] * 1 applies as long as 
e ~ 1 and d ~ 2/3 hold. 

However, the main dependence of a« on e and d is encased 
into b/{ = fan p v^/ksTn through the boundary temperature. We 



have just seen that T2 



v R A(p holds as long as strong 



shock conditions apply. If so, bR oc 1/A0 oc e would rise fast 
as e exceeds 1; but v 2 oc d 2 ^ 3 /e 2 is meanwhile reduced and the 
virial accretion shocks weakened, slowing down the net growth 
of b R . The overall result is that b R grows from the standard value 
2.7 to about 5 and then to about 8 as e increases from 1 to 1.5 
and then to 3. 

After Eq. (5), increasing values of b R cause or to decrease, 
and imply progressive saturation and even a decline of the en- 
tropy produced at the boundary. Specifically, as e increases from 
1 to 1.5 and then to 3, the entropy slope or goes from standard 
value 1 . 1 to about and then down to negative values around -2. 

To tackle the issue, recall that while the cluster outskirts de- 
velop to the currently observed radius R, the ICP is adiabatically 
compressed and settles into the DM potential well. Then the spe- 
cific entropy k(r) stratifies shell by shell, leading to a running 
slope a(r) = a R that re tains the sequence of the values set at th e 
time of deposition (see lTozzi & Normanl200 it lLapi et al.l2005l) . 
As a result, on moving out from the early cluster body to the out- 
skirts currently building up by secular accretion, the whole slope 
a(r) of the outer ramp decreases, and k(r) oc r fl(r) flattens or even 
bends over. 

To describe this behavior, in Lapi et all (1201 Oh we used an 
entropy slope a(r) = a — a' (r- n,) smoothly decreasing from the 
body value a ss 1.1. We found a value a & at r « R/2 an d val- 
ues a « -2 at R (as illustra ted in Fig. 1 of Lapi et all2010 ), con- 
sistent ly with the data by iBautz et aTl d2009t) and George et al.l 
(2009). 



2.3. Onset of turbulence 

Here we show that such an entropy demise due to M dwin- 
dling will arise together with the onset of turbulence triggered by 
shock weakening. In fact, the latter causes the postshock speeds 



V2 «i 1 3 1 
vi ~ «2 ~ 4 4 Al 2 ' 



(6) 



to grow from values around 1/4 and approach 1. So with V\ oc 
e~' lowering sharply while c\ varies as (1 + z) or less, the kinetic 
energy ratio is enhanced relative to the strong shock value 6.3%, 
see Sect. 2.1. For example, as e grows from 1 to 1.5 and then 
to 3 - 5, the ratio v 2 /v 2 increases from 10% to 14% and then 
to 25 - 39%. Note that the ratio v\/c\ of the residual bulk to 
the sound's speed past the shock also increases for decreasing 
M 2 ; in particular, it goes from 25% to 50% as M 2 ranges from 
1 to 3 - 5. On average, for a CC cluster we find the condition 
Ai 2 oc d 1 ^ 3 /(e(z))(l + z) < 3 for shock weakening to be met at 
redshifts z < 0.3, on using for the average (e) the values given by 
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Fig. 3. Profiles of pressure and temperature computed with the SM. The dashed line illustrates the pure thermal (laminar) case, while 
solid lines illustrate the turbulent case with S R = 40% and different t (color code as in Fig. 2). In the left panel, solid lines refer to 
the thermal pressure, while dotted ones refer to the total pressure. 



lLapi & Cavaherel {2009) in their Fig. 6. In a nutshell, combining 
Eqs. (4) and (6) yields an inverse relation of v^/Vj with kilk\. 

These postshock flows provide bounds to the energy level of 
subsonic turbulent motions that is driven by smooth accretion 
in relaxed clusters; minor, intermittent and localized contribu- 
tions may be added by the complementary clumpy component 
of the accretion, recently re-calibrated to less than 30% in the 

I II /J I I 

outski rts of relaxed halo s (see Fakhour i et al.H20 10: Gene l et al.l 
2010: [Wang et al.1 12010). Our bounds actually constitute fair es- 
timates for the amplitudes of outer turbulent energy at r ~ R, 
as shown by similar values obtained b oth observationally (see 
Mah davi et aDl2Q08t IZhang et al.l[20Toh and in man y numerical 
simulations from lEvrardld 19901) to lTau et al.l d2009l) . 



3. Modeling the turbulent support 

Turbulent motions start at the virial radius R with com- 
parable coherence lengths L ~ R/2, set in relaxed CC 
clusters by th e pressure scale height or by shock seg- 
mentation (see [lapichino & N iemeveri 120081, iRvu et alTI |2008, 
iPfrommer & Jonesl2010llVazza et al.l2010l) . Then they fragment 
downstream into a dispersive cascade over the 'inertial range', 
to sizes ( where diss i pation begins af ter t he classic picture 
bvlKolmogorovl (|1941|) . IObukhovl (ll941l) and iMonin & Yagloml 
(119651) . In the ICP context the dissipation scale (equivalently to 
the classic Reynolds' scaling) writes t ~ (c2/v) 3 ^ 4 A. pp (L/A pp ) 1 ^ 4 



in terms of the ion collisional mean free path A pp and of 
the ratio v/c? of the turbulen t rms speed to the sound's, 
see Unogamov & Sunvaevl d2003l) . For subsonic tu rbulence with 
vj_£2~ 1/3 (see direct observational bounds by ?, Sanders et al. 
1201 Ol and references therein) the relevant scale £ exceeds some- 
what A pp ~ 10 2 kpc. 

In the presence of outer magnetic fields B < I P" 1 uG (see 
Bona fedeet all 120101 IPfrommer & Jonesl 120101: also iRvu et al.l 
2008) the key if quantitatively debated feature is their degree of 



tangling; then then the effective mean free path is provided by the 

coherence scale of the field, somewhat larger than /l DD (see dis- 

I II 1 \^ ^ — * I 

cussions by Naravan & Medvedev 2001, Inogamov & Sunvaev 

l2003llGovoni et al.ll2006llBrunetti & Lazariarl2007l) . 

Clearly the above estimates provide useful guidelines, but 
need theoretical modeling and observational probing. Both is- 
sues are addressed next. 

Since turbulent mot ions contribute to the pressure (see 
lLandau & Lif shitzl [.1959) to sustains HE, our modeling is fo- 
cused on the ratio S(r) = p n th/Pth of turbulent to thermal pres- 
sure (or, equivalently, on the ratio 5/(1 + S) of turbulent to total 
pressure), with radial shape decaying on the scale I from the 
boundary value Sr. 

The total pressure is conveniently and generally written as 
Pth( r ) [1 + S(r)], while the thermal component is still expressed 
as pih(r) °c k(r) n 5 ^ 3 (r). With this addition, we proceed to solve 
the HE equation just along the steps leading to Eq. (1), and now 
find the temperature profile in the form 



T(r) 
Tr 



\k(r)} 


3/5 


1 + Sr 


2/5 | 


[ k R \ 




l+6(r) 





l + 2 - 



r R dx v 2 c (x) 

Jr X V 2 
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k(x) 



3/5 



5 1 +6 R 

3/5 



1 +6 R 



1 + 5(x) 



(7) 



which extends Eq. (1) for 5 > 0. Again, n(r) is linked to T(r) by 
n = [k B T/k] 3/2 . 

The actual temperature at the boundary is now lowered to 
Tr = T2HI + Sr). This is seen from Rankine-Hugoniot condi- 
tions in the presence of turbulent pressure; the latter obviously 
implies the term p2 (1 + Sr) to be added on the right hand side of 
the stress balance, and the corres ponding one 5 P2 (1 + Sr) V2/2 
in the energy flow (see Eqs. Bl in Lapi et al. 2005|). 
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In our numerical computations that follow we will adopt for 
fully developed turbulence the simple funct ional shape (rather 
than a constant 6 as used by Bode et al. 2003 

-(R-r) 2 /P 



6(r) = 5 R e 



(8) 



which decays on the scale ( inward of a round maximum, a 
smoothed out representation of the inertial range. This is illus- 
trated in Fig. 2 for three cases with £ - 0.9 R (blue), I = 0.5 R 
(red) and £ = 0.2 R (green). Note that the runs 6(r) we adopt 
are consiste nt with those rec ently indicated by numerical simu- 
lations (e.g jLauet al]|2009h . 

4. Results 

We illustrate in Fig. 3 the resulting temperature and pressure pro- 
files. Relative to pure thermal HE with 6r = 0, in the outskirts the 
thermal pressure lowers since it is helped by turbulent motions 
in sustaining the equilibrium, while the pressure run is moder- 
ately steeper; at the center the pressure is mainly contributed by 
the thermal component. 

As to the temperature (with the classic CC shape), it is seen 
that the variations are mild, and primarily stem from the reduc- 
tion of Tr/T2 by the factor 1 + Sr = 1.4 discussed in Sect. 3. 
Thus we recognize the saturation or bending over of entropy on 
scales r < R/2 (see Sect. 2) to constitute the primary cause for 
the steep temperature profiles observed by Suzaku; these drop 
by factors around 10 from the peak to the outer boundary in the 
low-z clusters like A1795 and PKS0745-191. 

This view is confirmed by Fig. 4, where we illustrate our 
bestfits with the SM to the temperature and surface brightness 
profiles for the two clusters A 1795 (north sector) and A 1689 
(azimutally averaged), one at a low and the other at a relatively 
high z. The figure shows that the SM with an entropy slope de- 
creasing through the outskirts from standard values a a 1.1 at 
r < R/2 to a k -1.8 at R (just as expected in Sect. 2.2) fits the 
temperature profile of the low-z cluster A 1795 much better than 
the case with uniform slope a * 1 . 1 . As expected, the difference 
is barely discernible for the relatively high-z cluster A1689. Both 
fits are only mildly affected on including in the SM the turbulent 
support. Note that the linked n(r) profiles flatten out to enhanced 
brightness landings (cf. insets in Fig. 4), a simple warning of 
interesting temperature and turbulence distributions. 

A far reaching consequence of turbulent support shows up 
(see Fig. 5, left) in the reconstructions of DM masses from 
X-ray o bservables bas ed on reversing the thermal HE equa- 
tion (cf. Sarazin 1988, pag. 92). It is seen that the mass re- 
constructed on ignoring the non-thermalized component devi- 
ates from the true mass by different amounts, and for small val- 
ues of the dissipation scale £ may even have a non-m onotonic 
behavior, not unlike the results by Kawaharada et al.1 (1201 Ol) in 
three sectors of A1689 (cf. their Fig. 8). In all cases the re- 
constructed mass within R is negatively biased by 20 - 30% 
relative to the true one. Accounting for such a bias consti- 
tutes a key point to resolv e the tension between weak lens- 
ing and X-ray masses (e.g., iNagai et aDl2007t lLau et al.| [2009; 
Meneghet ti et al.l 120 10). and in deriving precise cosmological 
parameters from st atistics of cluster m asses via the fast X-ray 
observations (see lVikhlinin et al. 2009). 

The above picture may be double-checked in individual clus- 
ters on directly gauging with the SZ scattering how the elec- 
tron pressures are lowered in the presence of outer turbulence. 
In Fig. 5 (right) it is seen that along l.o.s. running mainly into 
the outskirts we expect the latter to be lower by 20 - 30% rela- 
tive to the pure thermal case. This is stronger by a factor about 5 
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Abell 1689 @ z = 0.18, R 200 =2.9 Mpc 




in 0.0 0.2 0.4 0.6 0.8 1.0 



0.0 



0.2 



0.4 0.6 
S/R200 



0.8 



1.0 



Abell 1795 @ z = 0.06, R 200 = 1 .9 Mpc 




Fig. 4. Profiles of projected X-ray temperature (brightness in the 
insets) for the CC clusters A 1 689 ( top) and A 1795 (bottom). 
Data are from Snowden et al.l d2008l) with XMM-Newton (cir- 



Baut z et alj d2009l) and lkawaharada et al.l ((2010) 



cles), and from 

with Suzaku (squares). The solid lines represent our bestfits with 
the SM extended to include turbulence after Eqs. (7) and (8) 
with Sr = 40% and I = 0.5 R. The dashed ones illustrate the 
outcomes in the absence of turbulence, but still with entropy de- 
creasing outwards as in Lap i et alj d2010l) . For comparison, the 
dotted lines illustrate the case with uniform entropy slope. 



than the effects from delayed equipartition between electron and 
ion temperatures downstream the shock over the relevant mean 
free path. A corresponding reduction is implied in the require- 
ments for sensitivity and/or observation times with upcoming so - 
phisticated instruments like ALMA (see I Wong & Sara zin 2009). 
Meanwhile, statistical evide nce of SZ reductions has been ex- 
tracted from stacked data by Komatsu et al.l d2010h . 

5. Discussion and conclusions 

To answer the issue raised at the end of Sect. 1, we have in- 
vestigated the connection of ICP turbulence in cluster outskirts 
with virial accretion shocks; among the variety of shocks found 
in numerical simulations, these are best amenable to a simple 
treatment, that may also shed light on more complex conditions. 

A result to stress is the inverse nature of the connection that 
follows from combining the key Eqs. (4) and (6), namely, as the 
Mach number M of the shock decreases the postshock entropy 
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Fig. 5. Left: profile of DM mass; the dashed line illustrates the true mass, while the solid lines illustrate that reconstructed from 
X-ray observables on assuming pure thermal HE. Right: Projected profile of SZ scattering; the dashed line refer to the pure thermal 
case, while solid lines refer to the turbulent case. In both panels Sr - 40% and different values of I (color code as in Fig. 2) are 
adopted. 



k.2 oc Tijn^ lowers with Ai 2 while the residual bulk energy 
v\lv\ rises with M 2 . In other words, saturation of entropy pro- 
duction and increasing bulk flows to drive turbulence occur to- 
gether. 

To pinpoint when this is bound to occur, we hinge upon the 
quantity M v 2 {v2lv\) 2 oc d 5 ^ e~ 3 (V2M) 2 including three fac- 
tors: the infall speed v\ oc d 2 ^ /e 2 , see Eq. (3); the accretion rate 
M oc d/e; and the residual kinetic energy ratio v 2 /v 2 , see Eq. (6) 
and discussion thereafter. Such a quantity depends strongly on 
d and even more on e, which render the effects on the outskirts 
growth of the cosmological expansion accelerated by the Dark 
Energy, and of the declining shape of the initial DM perturba- 
tion wings, respectively. 

This leads us to predict that the turbulent support starts to in- 
crease on average for z < 1/2 during the late development of the 
cluster outskirts when d » 1/2 and e > 1 apply. Then turbulence 
briskly rises at z < 0.3 when the shocks become transonic. 

On the other hand, environmental (even anisotropic) vari- 
ance is introduced when the residual energy flow M v?, strongly 
dependent on e, is modulated by the adjacent filament/void struc- 
ture. At very low z a high level of thermal support will persist 
in sectors adjacent to filaments , as is apparently the case with 
A 179 5 and PKS0745-191 (see iBautz et aDl2009L iGeorge et all 
2009). By the same token, values e > 1 will prevail in cluster 
sectors facing a void, causing there an early onset of turbulence 
given sufficient M. Meanwhile, in neighboring sectors facing a 
filament values e w 1 still hold and the thermal support may 
still prevail, a condition that apparently ap plies to A 1689 (see 
iKawaharada et al.ll2010tlMolnar et al.ll2010h . 

With the closer focus provided by the Supermodel we find 
that as the shocks weaken not only the boundary entropy produc- 
tion saturates, but also the whole outer entropy distribution k(r) 
is to flatten and the temperature profiles T(r) to become steeper, 
consistently with the observations. Eventually, as the inflow ap- 



proaches the transonic regime, k(r) tends to bend over and T(r) 
to steepen somewhat further. 

But then subsonic turbulence arises at the boundary; its 
inner distribution follows the classic picture of the turbulent 
motion fragmentation, with a debated scale for its final dis- 
sipation. This demands closer probing, that we have tackled 
with the fast, analytic tool provided by the Supermodel. The 
present data from j oint X -ray and weak lensin g observations (see 
Kawaharada etafll201d: iMolnar et alj|2010h Izhang et alJ fcOlO) 



concur with simulations (see lLau et al. 2009) to show no evi 
dence for a dissipation scale much shorter than ~ 10 2 kpc, con- 
ceivably set by tangling of the magnetic field (see disc ussions by 
iNaravan & MedvedevfeOOlHBrunetti & Lazarianll2007l) . 

To summarize, in the outskirts of relaxed clusters we find: 

• Turbulence is related to weakening shocks at late cosmic 
epochs z < 0.3, when saturation of entropy production causes 
steep temperature profiles. 

• Turbulent excess pressures 6r = p nt h/ p± up to 40% arise at the 
boundary, declining inwards on scales t ~ 100 kpc. 

• The overall masses derived from X rays are necessarily biased 
low down to 20 - 30% when such a turbulent support is ignored. 

These findings are consistent with the current observational 
and numerical data. Moreover, we predict: 

• Variance concerning steep T(r) and turbulent support in cluster 
sectors will be correlated with the filament modulation of the 
adjacent environment. 

• The SZ scattering will be considerably lowered relative to the 
pure thermal case, along l.o.s. running mainly into the outskirts 
or their sectors where X rays concur with lensing data in signal- 
ing turbulent support. 

We stress that the SM provides a fast tool to represent and 
probe conditions of smooth inflows that prevail in the outskirts of 
relaxed CC clusters, away from mergers that scar the NCCs and 
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constitute the realm of detailed but time-consuming numerical 
simulations. 

A final comment concerns the giant radiohalos observed at 
the centers of several clusters to emit synchrotron radiation; 
these suggest that non-thermal support may be contributed by 
a mixture of magnetic field and relativistic par ticles acceler- 
ated by shocks and turbulenc e due to mergers (see l Brunetti et alJ 
I2007t iBiermann et all 120091: iBrunetti et all 120091) . These pro- 
cesses, with their limited acceleration efficiency and short per- 
sistence, are apparently widespread in NCC clusters, mainly at 
Z > 0.2; so they have minimal superposition or interference 
with the substantial, low-z, long-lived, outer turbulent compo- 
nent concerning mainly CCs, that we have addressed here. 

The picture we pursue envisages the infall kinetic energy to 
thermalize along two channels. First, supersonic inflows achieve 
thermalization at the shock transition via a sharp j ump (or a 
few j umps within a layer of limited thickness, see Lap i et alJ 
2005). Second, the subsonic turbulent motions left over down- 
stream the shocks, particularly downstream weak shocks, are 
dispe rsed into a cascade of m any effective degrees of freedom 
(see lLandau & LifshitzH l959). down to scales where dissipation 
becomes effective. The channels' branching ratio shifts toward 
the latter when the shock weaken, for z < 0.3; meanwhile, tur- 
bulence concurs with thermal pressure to support the ICP equi- 
librium in the outskirts. The picture substantiates the following 
formal remark: Eq. (7) corresponds to Eq. (1) for the variable 
T (1 + 6) in terms of the extended entropy k(l +6). 

This picture of onset and development of outer ICP turbu- 
lence may be fruitful in other contexts, in particular for shocks 
and turbulence in wakes around mergers. So it warrants close 
modeling with our fast yet precise analytic SM, to probe over 
wide cluster samples in a range of z the key turbulence features: 
amplitude and scale. 
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